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The Effect of a Single Supernova Explosion on the Cuspy 
Density Profile of a Small-Mass Dark Matter Halo 
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ABSTRACT 

Some observations of galaxies, and in particular dwarf galaxies, indicate a presence 
of cored density profiles in apparent contradiction with cusp profiles predicted by 
dark matter N-body simulations. We constructed an analytical model, using particle 
distribution functions (DFs), to show how a supernova (SN) explosion can transform 
a cusp density profile in a small-mass dark matter halo into a cored one. Considering 
the fact that a SN efficiently removes matter from the centre of the first haloes, we 
study the effect of mass removal through a SN perturbation in the DFs. We found 
that the transformation from a cusp into a cored profile is present even for changes as 
small as 0.5% of the total energy of the halo, that can be produced by the expulsion 
of matter caused by a single SN explosion. 
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1 INTRODUCTION 

The presence of a cusp in the centre of cold dark 
matter (CDM) haloes is one of the st r ongest results 
of N - body simulations ( Mo ore et alj 1 19991 : iNavarro et al.l 
1 19971 . Hool). However, the slope of this density profile 
is in apparent discrepancy with some observations of 
disk galaxies and g alaxy clusters, which exhibit rather 



flat density cores (Burkcrt 1995; Kuzio de Nara v et al 



20081 : Ide Blokl 120051: Ide Blok fc Bo sma 20021: Ide Blok et al 



2000) 



20031: iSalucci fc Burkertl 

Coccato et al.l (j2008j) studied the bulge and the disk kine- 



On the other hand, e.g.. 



matics of the giant low surface brightness galaxy ESO 323- 
G064 and showed that observations are not able to disentan- 
gle different density profiles, in that context. Parallel to the 
alternative models for dark matter already proposed, there 
are also attempts to solve this discrepanc y within the CDM 
cosmology u s ing baryonic physics ( e.g., iMashchenko et al] 
(|200rj . l2008l ): lGovernato et all (|2010l )). 

The so-called core/cusp problem has been in the spot- 
light of astroph ysical res e arch for quite a while (for a re- 
cent review, see Ide Blokl (|2010l )). One natural attempt to 
solve the apparent contradiction, is to study the gravita- 
tional effect of baryons in the dark matte r density pro- 
file. T his approach was first suggested by INavarro et al.l 
(1996). The study of external impulsive mass loss events 



jRead fc Gilmorell2005l ) and steady winds (jGnedin fc Zhaol 
2002) corroborated with the idea that baryons have a 
cruci al role in the det er mination of the d ensity pro- 
file dPedrosa et al J [20091 ) . IMashchenko et alj (|2006l ) and 
iPasetto et al.l (|2010l ) discuss how star fo rmation processes 
can in terfere in the central profile shape. IPenarrubia et ah! 
|2010l ) argue that the existence of the central cusp is nec- 
essary in order to maintain the baryons trapped in the 
halo's gravitational potential. According to their results, 
the sum of SN explosions, tidal effects and star formation 
processes can completely destroy dwarf spheroidal haloes 
which initially present a core like central density. Several 
authors suggested that the interstellar medium (ISM) of 
dwarf galaxy systems co uld be entirely rem o ved by super- 
novae (SNe) explosions jDekel fc Silkl Il98rj: iHensler et al.1 



novae (alNe) explosio ns (JJekel & bilk 1980; rlcnslcr et alj 
2004l:lMori et al.ll2002l.l2004ll997l : lMurakami fc Babull999l : 
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Mac Low fc Ferrara Il999l ). In this work, we study the effect 
of a single SN explosion on the removal of baryon gas from 
the first haloes and the effects of this removal on the shape 
of the density profile. 

As stated in IMashchenko et alj (|200rl l200Sfi once re- 
moved, the cusp cannot be reintroduced during subsequent 
mergers involved in hierarchical evolution of galaxies. This 
statement is supported by numerical and analytical results. 
Outco mes from simulations made by iKazantzidis et ahl 
(2006) imply that the universal characteristic shape of dark 
matter densit y profile s may be set early in the evolution 
of halos and IPehnenl l|2005h shows theoretically that the 
remnant cusp cannot be steeper than any of the progeni- 
tor cusps. 
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In this work, we consider the possibility to solve 
the apparent core-cusp prob lem via baryonic physics 
ijMashchenko et all [2009 . l2008h . The sudden gas removal 
makes the dark matter distribution to expand. Such mech- 
anism can operate efficiently in small high redshift haloes, 
creating dark matter cores (jCovernato et all l201Ch . Usin g 
arguments similar to thos e presented bv lTonini et al.l (|2006h ; 
Ide Souza fc Ishidal (|2010r ). we analyze this process in the 
context of distribution functions. 

It is initially assumed that baryons and dark matter 
are in an equilibrium configuration. Considering that energy 
perturbations track the initial energy distribution (domi- 
nated mainly by dark matter), after being perturbed, the 
halo will redistribute dark matter particles in a new equi- 
librium configuration. As a consequence, variations in the 
initial halo distribution function, caused by energy or angu- 
lar momentum perturbations, are followed by macroscopic 
quantities, such as the energy density profile. 

In this paper, we considered the ACDM model with best 
fit parameters from Ijarosik et all <|201ll ) ( WMAP- Yr^Q) , 
Q m = 0.267, fi A = 0.734, fl b = 0.0449, and H = 0.71. The 
paper is organized as follows. In section [21 we briefly review 
our knowledge about the propagation of blastwaves in the 
interstellar medium (ISM) and show the fraction of gas that 
can be removed from the haloes; in section [3] discuss the 
changes in density profile due to gas removal; and in section 
[4] present and discuss our results. Finally, in section [5] we 
discuss our conclusions. 



2 EXPULSION OF BARYONIC GAS BY A 
SUPERNOVA EXPLOSION 

Since the gravitational potential of the first collapsed 
haloes is shallow, the ionizing radiatio n from the first 
stars can expel the gas out of them (lAbel et al. 20071 : 



Alvarez et alt 12006); iKitavama et"afl |2004| ; IWhalen et al l 
2004 1 Wise fc Abe]||200sl ). As a result, a subsequent SN can 
break away from the halo due to the de creased gas density by 
photoionization prior to th e explosion l|Kitavama fc Yoshidal 
120051 ; IWhalen et all l2008af ). In order to understand the na- 
ture of the SN shock expanding into an essentially uniform 
and ionized intergalactic space, we assume spherical symme- 
try. 



2.1 Evolution of a supernova blastwave 

As described in the r eview on astrophysical blastwaves 
by lOstriker fc McKeel |T988), and in the recent paper of 
ISakuma fc Susal (|2009r i. the evolution of a supernova rem- 
nant (SNR) in the intergalactic medium (IGM) can be de- 
scribed as follows. Initially, the energy is largely thermal, but 
as the supernova expands, the adiabatic expansion converts 
thermal into kinetic energy. 

During the first stage, the SN ejecta sweeps out roughly 
the same amount of mass as its own in the surrounding 
medium. In the second stage, the expansion of the shock 
front is well approximated by the Sedov- Taylor solution. 



Eventually, radiative losses from the SNR become signifi- 
cant, and the remnant enters in the third, radiative, stage 
of its evolution. A thin shell is formed just behind the shock 
fron t. Finally the S NR ex pands conserving momentum. 

IVasiliev et al.l (|2008h found that supernova explosions 
with an energy 10 53 ergs expel a significant portion of the 
initial bar yonic mass from prot ogalaxies with total mass 
~ 10 7 Mm. lWhalen et all (|2008bl ) performed numerical sim- 
ulations of primordial supernovae in cosmological haloes 
from 6.9 x 10 5 - 1.2 x 10 7 M and showed that even less en- 
ergetic explosions are capable of ejection of more than 90% 
of the baryons from haloes containing < 1O 7 M0. Based on 
such studies, we assume that feedback from SN explosions 
are able to expel up to all the baryonic gas in primordial 
haloes. Our next step then, is to evaluate what fraction of 
the total halo mass is in form of gas. 



2.2 Gas mass fraction 

In order to evaluate the ratio of ejected gas mass to to- 
tal (virial) ha lo mass, we used the expression obtained by 
IGnedinl (|2000l ) 



f gaB (M,z) = 



1 + 0.26 



M 



(X) 



where f gas is th e mass fr action of the halo, Mf is th e 
filter mass (e.g . . IGnedinl l|2000h : iRodrigues et"aD (|201(t ): 
Ide Souza et al.i (|201ll )'). f2(, and Q, m are the baryonic and 
dark matter energy density parameter, respectively and M 
corres ponds to the tot a l (vir ial) mass of the halo. Accord- 
ing to iKravtsov et al.1 (|2004| ). M F k 5x 10 6 M at red- 
shift z = 10, and the reionization epoch is assumed to occur 
between z ~ 8 — 11 (e.g., the gas fraction of haloes with 
M ~ 10 7 M Q was f gas = 0.148 at z ~ 10). 

Since it is unlikely that the ejection of baryons was com- 
pletely efficient, we consider cases where 5%, 10% and 15% 
of the total mass was expelled, which correspond to 31%, 
61% and 92% of expulsion of the halo's baryonic mass, re- 
spectively. 



3 EVOLUTION OF THE DENSITY PROFILE 
OF DARK HALOES 

In order to understand the effect on the density profile, 
caused by the removal of baryonic matter from the halo, 
we study the evolution of its distribution function (hereafter 
DF). 

3.1 Distribution functions and density profiles 

The DF fully describes the state of any collisionless 
system at any time, specifying the number of particles 
f(x,v,t)d 3 xd 3 v having positions in the small volume d 3 x 
centered on x and velocities in the small range d 3 v centered 
on v. The evolution of a collisionless system is governed by 
the Boltzmann equation, ^7=0. 

The DF of a mass distribution in a steady state is re- 
lated to its density profile through 
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p(r) 



f(r,v)d v 



(2) 
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For an isotropic galaxy, the density can be written in 
the simple form 

p(r) = 4V2n [ * [*(r) - £} 1/2 f(£) d£ , (3) 



where ^ is the relative potential and 

£ = *(r) -v 2 /2 

is the relative energy (Cud defordlll99ll : iBinnev fc Tremaind 
l2008h . 

Equation ((3J is an Abel integral whose solution is 



4tt 2 d£ J y/£ - * 



(4) 



3.2 Parametrizing the density profile 

The usual parametrizations of cuspy density profiles (eg. 
NFW) lead to potentials without a simple analytical solu- 
tion for the inverse function, As a consequence, it is 
generally not possible to obtain an analytical form for the 
DF from equation Q. 

In order to avoid these complicat ions, we adopt the fo l- 
lowing family of spherical potentials (|Rindler-Da llcr 2009), 



* = 



(5) 



(fe Q + r a )T ' 

which leads to density profiles of the form 

, , .or, (1 + a)b a + (1 - cryK 

p(r) = Pgai b r2 , Q(ba+r , Q)7+2 , 

where p g is a characteristic density to be adjusted for each 
profile, or, in terms of the relative potential 



(6) 



1 + a V / 

x [l - aj+ a(l + 7)*^] 



This form of the density profile can conveniently capture 
the behavior of cusp and cored density profiles, for specific 
choices of the parameters (a, 7, 6). 



3.2.1 NFW profile 

The most common choice for the form of the density pro- 
file of haloes, which is also in best agreemen t with DM N- 
body simulations, is the NFW density profile (|Navarro et al.l 

E222), 



p(r) 



£ + 



(8) 



Where r s and p 3 are determined by the concentration pa- 
rameter c through 

3 

r s = and p s = — - — — , (9) 

c 3 log(l + c) - 

where p v [ T — 178 p, is the average density of a virialized halo, 
and r v i r corresponds to its radius. 

Since we are treating very high redshift (z > 10) haloes 
(i.e. in the beginning of their mass ac cretion histories), we 
use c = 4, following the prescription of IZhao et alj ((2009) . 
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Figure 1. In the upper panel, the continuous (gray) curve shows 
the NFW density profile of a halo of mass 10 7 M at z = 10, and 
the (blue) points correspond to the general profile with a = 1, 
7 = 1/2 and b = 1.085 r B = 36.8 pc. In the lower panel the 
fractional difference between the two profiles is shown. 



A different choice of c should not, however, affect sig- 
nificantly the forthcoming results, since all of them can be 
expressed in function of the parameters p s and r s and will, 
thus, scale with a change in the concentration. 

In figure [T] we show that the density profile of equation 
(0 provides a very good approximation to the NFW profile 
for the choice a = 1, 7 = |, b = 1.085 r s and p g — p s , with 
the difference between the two profiles never exceeding 10%. 

3.2.2 Cored profiles 

Many recent observations (|Burkertl Il995l ; iDonato et al.l 
2009) favor the Burkert density profile, which has the form 



Pbur(f) = 



(10) 



The Burkert profile can be well approximated by equa- 
tion ([6]) choosing a = 2, 7=^,6 = 0.76 ro and p g — po, as 
shown in figure [2] 

Another common choice for a cored density profile, is 
the pseudo-isothermal (PI) profile 

po 



Piso(r-) = 



1 + 



(11) 



3.3 Perturbing the distribution function 

The Boltzmann equation of a collisionless isotropic system 



d£ 
dt 



The removal of the baryons from the halo will lead to 
a small change in its energy, 8£. This change will take a 
small amount of time, 5t, and will cause a change in the 
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Figure 2. In the upper panel, the continuous (gray) curve shows 
the Burkert density profile of a halo of mass 10 7 M Q at 2 = 10, 
and the (blue) points correspond to the general profile with a = 
3/5, 7 = 6/5 and b = 1.58 r s = 53.7 pc. In the lower panel the 
fractional difference between the two profiles is shown. 



DF. We will assume that the distribution function, after the 
perturbation, f(t + 8t,£), will have the form 

f(t + 5t,£) = f{t,£ + 5£) . (13) 

To demonstrate the validity of the above assumption, 
we first expand the last equation keeping only first order 
terms in SE, 



f(t + 5t,£) = f(t,£ + S£): 



f(t y £) + 5£^ 



(14) 



If the ejection took place on a small time interval, St, the 
transformation in the distribution function can be written 



8f _ f(t + St,£)- f(t,£) _ S£df 
St St ' St d£ 



(15) 



which, in the limit of small time intervals, is precisely the 
Boltzmann equation. 



Relation between relative energy and 
perturbed halo mass 



3.4 



The formalism discussed in section 13.31 allows us to under- 
stand how the removal from the halo, of a certain amount of 
relative energy, S£, change the density profile. It is, however, 
necessary to relate this change of energy with the actual ex- 
pulsion of baryons by SNe. 

From a fundamental point of view, this correspondence 
is far from trivial since the same amount of energy could 
be removed by either lower mass, faster particles or higher 
mass, slower particles. There is, also, a cosmological context, 
where taking into account the mass outside the virial radius 
is not a sensible choice. 

We avoid those complexities taking the final density 
profile calculated for a given S£ and, then, calculating the 




r (pc 



Figure 3. New equilibrium configuration of the perturbed halo 
removing 5% (dashed curve), 10% (dotted curve) or 15% (thin 
curve) of the total mass of the halo, compared to the NFW profile 
(top solid line). 



virial mass by a simple 
procedure. 

We initialize the virial radius variable 



and fast convergent - iterative 

(0) 



with the 



virial radius of the unperturbed halo. We then calculate the 
virial mass using 



M 



.(i-l) 



47tr p(r) dr 



(16) 



where p(r) is the perturbed density profile. We reevaluate 
the virial radius 



M 



(i) 



4 

3 Pvir 



1/3 



(17) 



and proceed to the next iteration. 

This procedure converges rapidly to a consistent value 
for the virial mass of a perturbed halo, which can be com- 
pared with the known virial mass of the unperturbed halo 
(set to 10 7 Mq in our calculations) , therefore allowing us to 
relate variations in energy with variations in mass. 



4 RESULTS 

Using the formalism developed in the end of section [3l we 
calculated the DF associated with NFW density profile (the 
particular expression for the DF can be found in appendix 
[A| . The DF found were then perturbed through equation 
H13[) and a transformed density profile was generated from 
it using equation ([3]). 

In figure [3] we see how a NFW density profile varies 
due to the expulsion of 5%, 10% and 15% of mass of the 
halo (namely, removal of 0.57%, 1.17% and 1.81% of the 
relative energy, respectively). The appearance of a core in 
the transformed density profile is noticeable. 

In order to get better quantitative insight, we fit a 
pseudo-isothermal density profile to the resultant trans- 
formed profile. We found a core radii of ro = 0.054 r s , 
r = 0.075 r s and r = 0.092 r s for the removal of 5%, 10% 
and 15% of the the halo's mass, respectively. 

We also fit a Burkert-like profile finding ro = 0.20 r s , 
r = 0.24 r s and r = 0.28 r s for 5%, 10% and 15% of re- 
moval of the halo's mass. 
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Figure 4. New equilibrium configuration of the perturbed halo 
removing 1% (dotted line) and 2.5% (dashed line) of total energy 
E, compared to the original Burkert profile (solid line). 



but, instead, a core. We found that the transformation from 
a cusp into a cored profile is present even for changes as 
small as 0.5% of the total energy of the halo, that can be 
produced by the expulsion of matter caused by a single SN 
explosion. 
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As in the NFW case, we fit a pseudo-isothermal density 
profile to the resultant transformed profile. We found core 
radii of ro = 0.011 r s , ro = 0.012 r„ and ro = 0.013 r s for the 
removal of 5%, 10% and 15% of mass of the halo (namely, 
removal of 0.70%, 1.44% and 2.22% of the relative energy, 
respectively). The fit of a Burkert-like profiles leads to ro = 
0.074 r s , r = 0.069 r s and r = 0.064 r s for 5%, 10% and 
15% of removal of the the halo's mass. Ou r core radius are 
appro ximately similar to the one found bv lGnedin fe Zhiaol 
(2002). In their analysis, they ro lies in the range < ro < 
0.2r s , depending of the model for outflow. 

As the SN could also occur in a halo which was origi- 
nally cored (e.g. haloes where the core was previously erased 
by the process exposed), we performed the same analysis for 
a halo with an initial core profile. In figure f3] we show the 
evolution of a Burkert density profile after energy removal. 
As we can see, the core structure is maintained after gas 
removal and its radius remains approximately unchanged. 



5 CONCLUSIONS AND DISCUSSION 

We explored the effect of a supernova explosion on the cen- 
tral density of a small mass (~ 10 7 Mq) dark matter halo at 
high redshifts. 

We first reviewed (section [2]) the evidence that super- 
novae can efficiently expel a large part of the baryonic mass 
from small haloes. 

We then built a simple analytical model where we as- 
sume that 

(i) the halo is approximately isotropic, 

(ii) the mass expelled by the SN leads to a small loss of 
energy of the halo. 

Since typical parametrizations of the density profile do 
not lead to invertible expressions for the potentials, we used 
the profile of equation which is shown to be consistent 
with both cusp and cored density distributions. 

From the density profile we calculated the distribution 
function associated with NFW profile. We, then, evolved this 
distribution function, removing a small amount of energy 
from it. 

The final density profile found does not present a cusp, 



APPENDIX A: DISTRIBUTION FUNCTION 
NFW distribution function 

The DF found for the choice of parameters a = 1, 7 = | 
and b = 1.085r s , that emulates a NFW density profile, is 



(-1 + £ 2 ) 3 tv 2 



4(-l + £)£ 



3/2 



x (l+£)(28-9£ -136£ 4 96£: 
+ 2lV 1 — £ (1 + £ ) arcsin 



(Al) 



21(-1 + £)Vl + £ arctanh 



l + £ 
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